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Abstract 

How important would be a precise assessment of the electron screening effect, on 
determining the bare astrophysical S"- factor {Sb{E)) from experimental data? We 
compare the Sf,{E) obtained using different screening potentials, (1) in the adiabatic 
limit, (2) without screening corrections, and (3) larger than the adiabatic screening 
potential in the PP-chain reactions. We employ two kinds of fitting procedures: the 
first is by the conventional polynomial expression and the second includes explicitly 
the contribution of the nuclear interaction and based on a statistical model. 

Comparing bare 5-factors that are obtained by using different screening poten- 
tials, all Sh{E) are found to be in accord within the standard errors for most of 
reactions investigated, as long as the same fitting procedure is employed. Sh{E) is, 
practically, insensitive to the magnitude of the screening potential. 

Key words: Astrophysical 5-factor; Electron screening in the laboratories 



1 Introduction 



On the main- sequence stellar sites the series of reactions that convert hydrogen 
into hehum is known as the proton-proton chains. It is a key to understand 
the evolution of the stars. The cross sections of these charged particle induced 
reactions are major ingredients to calculate thermonuclear reaction rates. They 
are measured at laboratory energies and are then extrapolated to thermal 
energies (jH), because of their smallness at such low energies. The extrapolation 
is performed by introducing the astrophysical S'-factor: 

S{E) = a{E)Ee^''''^^\ (1) 
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where (t{E) is the reaction cross section at the incident center-of-mass (cm.) 
energy E and ri{E) = ZxZpaJ^, with Zt, Zp, and yU denotingthe atomic 
numbers and the reduced mass of the target and the projectile (y); a and c 
are the fine-structure constant and the speed of hght, respectively. The ex- 
ponential term in the equation represents the inverse of the Coulomb barrier 
penetrability. Since we have factored out the strong energy dependence of cr{E) 
due to the barrier penetrability, the S'-factor can be approximated by a smooth 
polynomial expansion in the absence of low-energy resonances. In laboratory 
experiments, the targets are usually in gas or solid state. In the low-energy re- 
gion, the S'-factors obtained from experiments show large enhancement to the 
extrapolation from high energy data for various reactions (jsl). This enhance- 
ment is, usually, attributed to the screening by the bound electrons around the 
target. In contrast in the stellar nucleosynthesis nuclei are almost fully ionized 
and surrounded by the plasma electrons. In deuteron induced reactions on 
deuterated metals (0; S M) and proton-induced reactions on lithium isotopes 
in several forms of lithium chemical compounds (0) , much larger screening en- 
hancements have been observed with respect to the enhancement by gaseous 
targets. The nuclear reactions in such a circumstance are affected by a differ- 
ent mechanism of the plasma or the conduction electron screening (@; [ol; Eoh . 
A similar effect has been discussed in the radioactive decay of a nucleus in a 
model (jlll ). The screening effects in the medium can be dependent on temper- 
ature and density of the medium and we do not consider such effects in this 
paper. Hence the screening effect of the bound electrons should be removed 
from the S'-factor data to asses the reaction rate in the stellar site correctly. 
The enhancement by the bound electrons is discussed in terms of a constant 
potential shift (screening potential Ue). The upper limit of Ue is obtained, 
when the adiabatic approximation is fulfilled, and it is given by the difference 
of the binding energies of the target atom and the united atom. (jl2l ) On this is- 
sue, dynamical effects have been studied by followin g th e time evolution of the 
atomic wave function in the classical allowed region (jl3l ). They solved the time 
dependent Hartree-Fock equation and evaluated the screening potential. Their 
results suggest that the screening potential approaches the adiabatic limit as 
the incident energy becomes lower. The influence of the tunneling phenomenon 
to this problem has been studied as well (jl4l). And, there, the screening po- 
tential could go over the above-mentioned adiabatic limit slightly, only in the 
case where the electronic wave-function has some excited state components 
at the classical turning point of the inter-nuclear motion. We have examined 
the problem using molecular dynamics approach with constraints (llSl: Il6l). to 
see the effect of the fluctuations ally . llSl ). The obtained average enhancement 
factors do not exceed the adiabatic limit. However, there are events that give 
enhancement factors larger than that in the adiabatic limit. 



In this paper, we discuss the influence of the electron screening effects on the 
determination of the bare S'-factors {Sb{E)). We especially focus on the reac- 
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tions in the hydrogen burning process and determine Sb{E) through a fitting 
of experimental data making use of the polynomial expression and fixing the 
screening enhancement in the adiabatic limit. There are several more sophis- 
ticated theoretical models that describe the bare S'-factor as the i?-matrix 
model (19; 2^; 21), the potential model (22), and the distorted wave Born ap- 
proximation. However our aim is to clarify the effect of the electron shielding 
on the experimental data of the ^-factors rather than to determine the S'-factor 
based on a particular theoretical nuclear model. For this purpose we choose 
the simplest way: the polynomial expression and try to determine Sb{E) in 
a model independent way. We will discuss the sensitivity of the 5*^(0) values 
on the choice of the degree of the polynomial and the fitting range for each 
particular reaction. We stress that the use of simple polynomial expressions 
does not carry physical meaning and hence the derived Sb{0) values should be 
taken with caution. To make up this deficit, we propose another expression in 
place of the polynomial expression. The new expression includes an explicit 
contribution of the nuclear interaction. Moreover it is based on a two-step 
process with a compound nucleus (CN) formation and a statistical choice of 
the exit channel (1231 : |24J : l25l : |26| ). Provided that the energy regions of the as- 
trophysical interest is low, we assume that the / = partial-wave component 
is dominant for most of the reactions. To our knowledge, the statistical model 
calculation has been used to estimate the S'-factor of the radiative proton- 
capture reactions on Sr isotopes (1271 ). They have used the Hauser-Feshbach 
statistical model code and have compared the theoretical results with the ex- 
perimental data. Because in the Hauser-Feshbach statistical model code, a 
global level density, which is based on some models, is used, one can obtain 
the absolute value of the S'-factor. They have found discrepancies between the 
theoretical results and the experimental data, especially in the reactions with 
target isotopes near the neutron shell- closure. In their paper they did not aim 
to fit their experimental data. 



Electron screening effects have been studied on some specific reactions in the 
chain, especially on transfer reactions ^He(^He,2p)^He and ''Li(p, a)^He, which 
are studied including extremely low-energy region (1281 : l29l : l30l ) . However there 
is no systematical study for all of the reactions including the radiative capture 
reactions. We aim to deduce these effects from such well studied reactions and 
estimate the screening quantitatively on all the reactions in the chain. This 
motivates us to employ the adiabatic limit, rather than leaving the screening 
potential as a fitting parameter as it is customary treated in a series of stud- 
ies ( 1281 : 120I ) . For comparison we also perform the fitting procedures without 
enhancement factor and by treating the screening potential as a parameter. 
The obtained Sb at zero incident energy are compared with the results in the 
NACRE compilation JH), in which the authors employed the screening po- 
tential higher than the adiabatic limit, and with the results in the i?-matrix 
analyses (12 ll ). 



3 



This paper is organized as follows, In Sec. [2] we describe the enhancement 
factor by the bound electrons within the adiabatic limit briefly. We list up the 
reactions in PP-chains in Sec. [3] and explain how we incorporate the enhance- 
ment factor into the fitting procedure. We, especially, give a detailed account 
on the second fitting procedure based on the statistical model. Some reactions 
are analyzed in this section. We summarize the paper in Sec. |H 



2 Enhancement factor in the adiabatic limit 



To discuss the enhancement quantitatively, we determine the enhancement 
factor (nj 



in terms of the measured cross section (j{E) and the bare cross section ao{E). 
If one assumes that the effect of the electron screening can be represented 
by the constant shift t/g (screening potential) of the potential barrier, the en- 
hancement factor is approximated by (jl2l : Il3l ) , 



E 

7n]{E) 



log fe- 



(3) 



The Ue can be estimated easily in two limiting cases. In one case the inter- 
nuclear velocity is much higher than that of electrons velocity; this limit is 
called the sudden limit. Within this limit the electron wave function is frozen 
during the reaction. In the opposite case where the inter-nuclear motion is 
much slower than electrons motion, the bound electrons follow the motion of 
the nuclei adiabatically. Within this adiabatic limit the screening potential 
is expressed by the difference of the binding energies between the initial tar- 
get atom (BEt), and the united atom (BEua), which is formed during the 
reaction. 

UiAD) ^ _ ^^^^ 



The screening potential within this limit gives the theoretical upper limit. 
However, we should stress that these cases deal with the ideal situation of 
a ion impinging on an isolated atom. This should be the ideal situation in 
nuclear physics experiments where very thin targets are used together with a 
well collimated mono-energetic beam. However, this situation is not fulfilled 
if the beam energy is very low (which is the case of interest in the present 
paper) or if the atoms are embedded in a medium such as a metal at a given 



4 



Table 1 

Reactions in PP-chains, their minimum incident energy in the cm. system measured 
so far and the enhancement factor within the adiabatic hmit at the minimum energy. 



reactions 




(keV) 


je \J^min ) 










D(p,7)3He 




2.52 


1.07 


3He(3He,2p)^He 




20.76 


1.22 


=^He(a,7)7Be 


93.#, 


127.* 


1.02 


^Be(e-,z^e)'^Li 








^Li(p, a)^Re 


12.7 


10.** 


1.18 


^Be(p, -ffB 




115.6 


1.01 



*prompt-7 method *activation method **THM 

density and temperature. In the latter case the nuclear process is expected to 
be influenced by the rearrangement of the electrons in the metal which will 
give rise to some peculiarities similar to the studied radioactive decay in a 
medium (llll). 



3 Bare S-factors of PP-chain reactions 



A list of reactions in PP-chains is shown in Tabled! The first reaction H(p, /3"^i^e) 
involves the /?-decay and has too small cross section to be measured experi- 
mentally. Its S'-factor is calculated from first principles (1311). We, therefore, 
concentrate on the other 5 reactions except the electron capture reaction 
''Be(e~,z/e) ''Li. In the table the minimum incident energies, measured so-far, 
for each reaction are also shown. For two transfer reactions, ^He(^He,2p)^He 
and ^Li(|}, a)'^He, the cross sections have been measured already including the 
low-energy region where the screening enhancement becomes more than 10%. 
The other three reactions are radiative-capture reactions, which have even 
smaller cross sections. The S'-factor of the reaction ^He(a,7)^Be has been re- 
determined with high precision recently both by detecting 7-ray from ^Be 
decay(the activation method) (1321 : l33l ) and by detecting prompt 7-ray(the 
prompt method) (jj^ ). Its S-factor in the low-energy region is extrapolated 
from high energy data by the i?-matrix fitting. The reaction ^Be(p, 7)®B in- 
volves unstable nuclei. The S'-factor of this reaction has been determined by 
means of the direct capture reaction (1351 : l36l ) and the Coulomb dissociation 
method (1371 ). It has been claimed that there is an inconsistency between the 
results of the two methods (i36, : i38i) . though the question seems to be resolved 
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by reanalyzing the data of the Coulomb dissociation method (1391 ) . 

If one takes into account modifications due to the nuclear potential and all the 
contributions from partial-waves, the bare S'-factor can be expressed as (0; |26l) 



Sb{E) 



2fi 



(5) 



where Ilif{E) is the probability to obtain a specified exit channel / from a 
certain entrance channel I and 
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AZrZpe'^ 
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2E 



sm 
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with E, = ZTZpe^R (MeV), Ei = l{l + l)h^/{2^iR^) (MeV fm^), i.e., the 
heights of the Coulomb barrier and the centrifugal potential at the nuclear 
interaction radius R. We assume an empirical formula R = d x Rnq, where 

1 /3 1 /3 

i^ATg = 1.4 X {Aj, + Ap ) (fm) and ci is a parameter that takes into account 
the fact that the nuclear potential has a diffuseness; with At and Ap denot- 
ing the mass numbers of the target and the projectile nuclei. We call d radial 
parameter. The exponential term in Eq. ([5]) stands for the penetration fac- 
tor divided by the pure coulomb penetrability. At zero incident energy limit, 
Eq. (jSD reduces to (1261 ) : 



2/i 



2l(l+l)H 



no/(0) + 5:%(0)(2/ + l)e 



i>i 



•(7) 



The conventional polynomial expression for the bare S'-factor of non-resonant 
reactions is associated to the Taylor expansion of Eq. ([5]). Most reactions in 
the PP-chain are non- resonant, one, therefore, can fit the bare S-factor using 
the enhancement factor: 



S{E) = Sh{E) ■ /e; StiE) = 3^(0) + S,E + S2E'' + ■ ■ •, (8) 

in an implementation of the nonlinear least-squares algorithm. For resonant 
reactions the energy dependence of the S'-factor is given by the Breit-Wigner 
formula (0) 

S{E) = — ^^^^^^^^ e^'^^(^) (9) 
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where uo = (2J + l)/(2ji + l)(2j2 + 1) and Er are the statistical factor and the 
resonance energy, respectively; ri(£'),r2 and F are the entrance and the exit 
channels partial widths and the total width. The incident energy dependence 
of Ti{E) is given again by the penetration factor in the case of sub-barrier 
reactions. One, therefore, can write down 



nh'^ Cre^^Uif{E) 
2/i {E - Erf + (r/2)^ 



where = coT2^^J^^ (MeV^), but we determine from the fitting proce- 
dure. In the case where there are more than one data sets we weight the data 
depending on its standard error. As one can easily imagine, higher order terms 
are important to fit the experimental data in the high incident energy region. 
In this paper we limit the fitting energy range less than 1 MeV and choose 
the degree of the polynomial to obtain convergence of the ^^(O) value within 
the statistical errors. We will discuss the sensitivity of 5*6 (0) on the degree of 
the polynomial and the fitting range for each particular reaction. 

Alternatively, the experimental data are fitted using Eq. ((51) directly, instead 
of the polynomial expression. For this purpose we derive the incident-energy 



dependence of Ilif{E) in Eq. ([5]). According to Weisskopf model (1231 : |24 l26l ). 
the probability to obtain a certain exit channel / after the CN formation is 
proportional to 

TTj oc gj{T',^ + 2mfn^) exp (-^=) <^ (H) 



where gf, Tk^, and mf are the number of states for the spin, the kinetic energy, 
and the mass, respectively, of the lightest reaction product in the exit channel; 
QcN; o'j''*, and a are the Q- value for the CN formation, the cross section of the 
inverse process, and the level density parameter a ~ Acn/8.0 (MeV~^)(2^; 
with AcN denoting the mass number of the CN. The absolute value of 11// (i?) 
is, in principle, given by 

^ifiE) = (12) 



where the sum in the denominator is taken over all possible exit channels. 
We, however, avoid calculating the sum and determine n;/(0) from fitting 
procedures of the experimental data. Instead, we scale the incident energy 
dependence of Ilif{E) 



^ _n ,^, n^{Ef + 2m,n^{E) (nAow^ ti^^WoA 



,(13) 
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where Tk^{E) = '^^^ ^ (E+Q), with Q, and A/ denoting the reaction Q- value, 
and the mass number of the hghtest reaction product. Eq. ([6]) is obviously valid 
only at the incident energy lower than the Coulomb barrier. In the fitting 
procedure using Eq. ((51), n//(0) (/ = 0, 1, 2, • • •) and the radial parameter d are 
treated as fitting parameters. We assume that the nuclear interaction radius 
can differ for each partial- wave: Ri = di x R^q. The fitting procedures are 
performed only in the energy region below the barrier. The sum of the partial- 
waves is taken up to the order with which the fit converges, mostly / = 0, 1, 2 
are sufficient. We anticipate that the description of the reaction mechanism 
through the CN formation works well, in particular, in the reactions with the 
/ =0 partial-wave in the entrance channel. (l26l ) 



3.1 3He(3He,2p)^He 



The ^-factor of the reaction '^He(^He,2p)^He from several measurements are 
shown with error bars in Fig. [H At the minimum incident energy, which has 



been reached in an experiment by the LUNA collaboration (|28| ). the screening 
enhancement is estimated to be more than 20% of the adiabatic approxima- 
tion. In the NACRE ([]:) compilation, the authors used the screening potential 
Ue = 330 (eV) and a quadratic polynomial to obtain the S'-factor. The fitting 
parameters in ([ll) are shown in the first row of table 2. 



For the reaction ^He(^He,2p)^He, the adiabatic screening potential is obtained 
under the following considerations. In the target medium ^He projectiles are 
likely to be ^He"*" or ^He charge neutral state. If we consider ^He neutral 
projectiles, the adiabatic screening potential is t/^^^) = 246.8 (eV) 0). 

rj(AD) 

/i^^)(=^i/e)=e-^(^)^. (14) 



For ^He"*" projectiles the adiabatic screening potential is calculated taking into 
account the charge symmetry of the system Uj^^^^^ = 255.5 (eV) and 

^{AD)2 ^ ^22.2 (eV) in the cases where the system ends up with ^Be"'"(ls)^(2s) 
state and ^Be"'"(ls)(2p)^ state respectively. The corresponding enhancement 
factor within the adiabatic limit is written as f.14) 



Txr]{E) 



E 



+ exp 



TT 



JJ{AD)2 

E 



(15) 



The results of the polynomial fitting using a quadratic polynomial are shown in 
Table [2] together with fitting parameters in yj). The obtained fitting parame- 
ters from the fit without enhancement factor are shown in the 4th row. The pa- 
rameters in the second row are for ^He neutral projectile and ones in the third 
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0.01 



0.1 
E[MeV] 



Fig 
the 



1. The S-factor for the reaction ^He(^He,2p)'*He as a function of 



incident CM. energy. The experimental points are from (j4Cl) (Backer67) , 
from ( 4ll) (DwarakanathTl ) . from (j42l )(Krauss87) and from ( 28l) ( Junker98) . The 
sohd curves represent the S'-factors multipUed by the enhancement factor within 
the adiabatic Umit from our fitting. The dashed curves show the corresponding bare 
S-factor. In the top and bottom panels the fitting is performed using the polynomial 
expression and using Eq.([5]), respectively. 



row are for ^He^ projectile. The corresponding curve for ^He neutral projectile 
case is shown in the top panel of Fig. [1] together with the experimental points. 
Notice that the adiabatic limit gives a smaller Sh{0) but within the standard 
error of the one obtained by the NACRE collaboration. If we fit the same 
data by varying the screening potential Ug, we obtain Ue = 371 ±46 (eV) and 
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Table 2 

Fitting parameters of the reaction ^He(^He,2p)^He using polynomial expression. 
The first row is from (Ij). The second and the third rows are obtained by using 
the screening potential in the adiabatic limit. The last row is obtained by assuming 
without enhancement. 



5b(0)(MeVb) Si{h) 52(MeV-ib) C/e(eV) xl 





5.18 






-2.22 




0.80 




330 




/i^^)(^He) 


5.23 


± 


0.06 


-3.1 ± 


0.5 


1.6 ± 


0.5 


246.8 


0.7 


/i^^)(3He+) 


5.32 




0.06 


-3.5 ± 


0.5 


1.9 ± 


0.6 




0.8 




5.56 




0.07 


-4.7 ± 


0.6 


3.0 ± 


0.7 





1.1 



Table 3 

Fitting parameters and the Sb{0) of the reaction ^He(^He,2p)^He using Eq. ([5]). 
The first row is obtained by using the screening potential in the adiabatic limit, the 
second row without enhancement, and the last row is obtained by treating Ue as a 
fitting parameter. 



no/(o) 


ni/(o) 


do 


di 


UeieY) 




5b(0)(MeVb) 


0.19 


7.4x10"^ 


0.61 


9.9 


246.8 


0.68 


5.4 


0.17 


2.2x10"*^ 


0.63 


13. 





0.73 


6.5 


0.18 


1.0x10-"^ 


0.60 


9.1 


299±102 


0.68 


5.2 



Sb{0) = 5.06±0.09 (MeVb) , with xl = 0.7. The deduced is higher than the 
adiabatic limit but the obtained S'b(O) is very close to the one in the adiabatic 
limit. Fixing the fitting range from the lowest experimental data 0.0208 MeV 
to 1 MeV, we obtained 5'f,(0)=5.32 ±0.08 (MeVb) using a cubic polynomial 
with x^=0.7. This ^^(O) coincides with the result using a quadratic polyno- 
mial (the second row in Tab. [2]). Limiting the fitting range from 0.0208 (MeV) 
to 1 (MeV), ^^(O) is insensitive to the choice of the degree of the polynomial. 
Assuming the quadratic polynomial and the adiabatic enhancement factor, the 
value SbiO) varies from 5.23±0.06 (MeVb) to 5.26±0.3 (MeVb), as one changes 
the upper-limit of the fitting from 1 (MeV) to 0.1 (MeV) but fixing the lower 
limit 0.0208 (MeV). On the other hand Sb{0) varies from 5.23 (MeVb) to 
5.17 (MeVb), as one changes the lower limit from 0.02 (MeV) to 0.2 (MeV). 
Thus the Sb{0) obtained by using the polynomial expression is not much sen- 
sitive to the choice of both the lower and the upper limit. 

The fitting procedures using Eq. ([5]) with the screening potentials Ue = 246.8, 
(eV) and treating Ue as a fitting parameter give zero energy S'-factors 5.4, 
6.5, 5.2 (MeVb), respectively, and they are shown in Table [31 We have used 
the / =0 and 1 partial-waves. The fitting procedure with Ue = gives xl 
slightly larger than the other two cases. If we take into account the enhance- 
ment factor, the radial parameters do and di for two cases are about 0.6 and 
9., respectively, and nij(O) is much smaller than no/(0). The former implies 
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Table 4 

Fitting parameters of the reaction ''Li(p, a)^He. The three rows are obtained by 
using a cubic polynomial and from ([l|) (the first row), with the adiabatic approxi- 
mation(the second row), and without enhancement (the third row). 



5^(0) (MeVb) 


Si{h) 


52(MeV-^b) 


53(MeV-2b) 






0.0593 


0.193 


-0.355 


0.236 


300 




0.0620 ± 0.0006 


0.15 ± 0.01 


-0.24 ± 0.03 


0.14 ± 0.03 


175 


0.41 


0.0673 ± 0.0008 


0.10 ± 0.01 


-0.13 ± 0.04 


0.08 ± 0.03 





0.62 



Table 5 

Fitting parameters and the S'b(O) of the reaction '^Li(p, Q!)^He using Eq. ([5]). The 
first row is obtained by using the screening potential 300 eV, the second row in the 
adiabatic limit, the third row without enhancement, and the last row is obtained 
by treating Ue as a fitting parameter. 



ni/(o) 




di 


Ue (eV) 




S;,(0)(MeVb) 


7.1x10- 


-5 


4.7 


300 


1.4 


0.033 


6.0x10" 


-5 


4.9 


175 


2.1 


0.035 


5.0x10- 


-5 


5.1 





3.5 


0.038 


9.6x10" 


-5 


4.3 


495 ± 41 


1.0 


0.031 



that the effective radius of the / =1 partial- wave component is larger than that 
of the / = component. And the latter implies that the / = component gives 
a dominant contribution to the S'-factor. The / =1 component plays a major 
role in the higher energy region. The resulting Sf,{0) are in agreement with 
the extrapolations using quadratic polynomials for the corresponding screen- 
ing potentials. The curve obtained by this fitting procedure for the adiabatic 
enhancement is shown in the bottom panel in Fig. [2l If we fit the same data 
by varying the screening potential Ue, we obtain Ue = 299 ± 102 (eV) and 
5*6(0) =5.2 (MeVb), with xl = 0.68. This xl is the same as in the case where 
we used the adiabatic screening potential. 



3.2 ^Li(p,a)^He 

In Fig. [2] experimental data of the S'-factor of the reaction '^Li(p, a)^He from 
several direct measurements are shown with error bars. We performed the fit of 
the data in the incident energy region from 0.01 MeV to 1 MeV using a cubic 
polynomial without enhancement factor. The obtained fitting parameters are 
shown in the third row of Table HI This fitting procedure without enhancement 
factor is quite sensitive to the choice of both the upper and the lower limits 
of the fitting range. 
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E[MeV] 



Fig. 2. 5'-factor for the reaction ^Li(p, a)'*He as a function of the incident cm. 
energy. In the top panel experi mental points are taken from (H) (Cassagnou62) . 



jxperi 

from (j45|)(Rolfs86), from ([291; l30|)(Engstler92). The sohd curve represents our fit by 
polynomial expression with the adiabatic enhancement. The dashed curve corre- 
sponds to the bare 5'-factor. In the bottom panel the data from (|4^)(Lattuada01) 
are shown, together with other experimental data in the top panel. The solid curve 
represents our fit by using Eq.Q, instead of the polynomial expression, and with 
treating the screening potential parameter. 
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The experimental data from (j29|; l30l ) show the enhancement of the S'-factor 
in the low-energy region. In the compilation NACRE ([l|), authors used the 
screening potential f/g = 300 (eV), which is larger than the adiabatic limit 



and a cubic polynomial to obtain the S'-factor. The fitting parameters in yj 



are shown in the first row. However in the experiments (1291 : ISOl ) LiF solid tar- 
gets and deuteron projectiles as well as deuterium molecular gas targets and Li 
projectiles are utilized. In the case of LiF target, which is an ionic crystal simi- 
lar to NaCl and a large band gap insulator, one can approximate the electronic 
structure of the target ''Li state by the ^Li"*" with only two innermost electrons. 
Thus one expects the screening potential in the adiabatic limit f/J^^^ = 371.8- 
198.2^174 (eV). The best fit of the experimental data in the form of Eq. 



and t/i^^) = 175 (eV) disl) gives the fitting parameters in the second row of 
Table m with the reduced = 0.46. The corresponding S'-factor is shown with 
the dashed curve in Fig. [2J If we fit the same data by varying the screening 
potential f/e, we obtain Ue = 195 ± 28 (eV) and ^^(0) =0.061±0.001 (MeVb), 
with xl = 0.41. This is, practically, consistent with the result (the second 
row of Table H]) in the adiabatic limit. S;,(0) in the adiabatic limit is lower 
than the value obtained by the fit assuming no screening enhancement but 
higher than the value obtained in (jl|) where the authors used Ue = 300 (eV). 
We have checked the sensitivity of S'b(O) on the degree of the polynomial 
and the fitting range. Fixing the fitting range from the lowest experimen- 
tal data O.OllMeV to IMeV, we obtained 5b(0)=0.0620±0.0006 (MeVb) and 
0.0617±0.0009 (MeVb), using a cubic polynomial and using a quartic polyno- 
mial, respectively. Both have x^=0.41. Hence we choose a cubic polynomial 
for the following fit. In addition to this, assuming the adiabatic enhancement 
factor, the value S;,(0) varies from 0.0620 (MeVb) to 0.0632 (MeVb), as one 
changes the fitting upper limit from 1 (MeV) to 0.5 (MeV) with fixing the 
lower limit 0.011 MeV. On the other hand S'b(O) varies from 0.0620 (MeVb) to 
0.578 (MeVb), as one changes the lower limit from 0.01 (MeV) to 0.04 (MeV). 
S'b(O) obtained by using the polynomial expression is more sensitive to the 
choice of the lower limit than to the choice of the upper limit. 

The fitting procedures of the same data but using Eq. are performed. 
We have used only the / =1 partial- wave, because / must be odd to obtain 



positive-parity state of ^Be ( 1471 : l2ll ). Using only the / =1 component, our 
fitting procedure gives a steeper incident energy dependence than the case 
where / =0 is used. In passing we mention briefly that the fitting procedure 
of the reaction ^Li(p,a)'^He in the next subsection for a comparison. With 
the screening potential Ue =300, 175, (eV) we obtain zero energy S-factors 
0.033, 0.035, 0.038 (MeVb), respectively, as they are shown in Table O They 
are all considerably smaller than the results of polynomial fitting and xt '^^ 
cases are larger than the results of polynomial fitting. The radial parameters 
di for all three cases are considerably larger than 1. This fact suggests that 
the interaction radius is 5 times larger than the empirical formula for the / =1 
partial-wave. Again, if we fit the same data by varying the screening potential 
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Ue, we obtain Ue = 495 ±41 (eV). The fitting parameters of tliis procedure are 
shown in the last row in Table [51 The curve obtained by this fitting procedure 
is shown in the bottom panel in Fig. [21 The extracted bare S'-factor data by 



THM are especially shown with the closed squares (l46l ) but these data are 
not included in the fitting procedure. Nevertheless the obtained bare S'-factor 
curve follows the data by THM, which is thought to give the bare S-factor. 

In Ref. (i2lt ) the i?-matrix fitting for higher energy region (E > 40 keV) has 
been used to determine the S'-factor. They obtained the zero-energy S'-factor 
0.067 ± 0.004 (MeVb) and the screening potential Ue = 100 ± 25 (eV), which 
is less than that within the adiabatic limit. The S'-factor at zero energy from 
our results using polynomials 5f,(0)=0.065 ± 0.005 (MeVb) and 0.0620 ± 
0.0006 (MeVb) are in agreement with this result from the i?-matrix fitting. 



3.3 6Li(p,a)3He 



In contrast to the reaction '^Li(p,a)^He, the reaction ^Li(p,a)^He does not 
have the restriction on the incident partial-wave. For a comparison, we show 
the fitting curves of this reaction in Fig. [31 where the top panel shows the two 
curves from a fitting procedure using a polynomial and the bottom panel shows 
the curves obtained by using Eq. ([5|), although this reaction is not included in 
the PP-chains, 

The fitting procedures of the experimental data using Eq. ([5|) are performed. 
Using the / =0 partial-wave alone, we obtain Sfe(O) = 3.3 (MeV), with xl =2.5 
with fixing the screening potential Ue = 175 (eV) in the adiabatic limit. We 
obtain Ue = 466 ± 31 (eV) and SbiO) = 3.1 (MeV), with xl =1-4 by treating 
Ue as a fitting parameter. The fitting parameters are shown in Table [HI For this 
reaction the radial parameter do for all three cases are about unity. Although 
the screening potentials are different, the difference does not affect either the 
Sb{0) or the fitting parameters, no/(0) and do. Only the xl ^^i Ue = 466 
±31 (eV) is much smaller than the others. The obtained screening potential 
from the latter procedure is in agreement with one in the reaction ^Li(p,a)'^He. 
This fact supports the isotopic independence of the electron screening. 

From the tree results of all the transfer reactions considered, ^He(^He,2p)^He, 
^Li(p, a)^He, and ^Li{p,a)^B.e, one can say that the enhancement by the 
screening is crucial, in the sense that the fitting procedure without enhance- 
ment gives xl larger than the others. However the obtained Sb{0) is insensitive 
to the magnitude of the screening potential. 
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Fig. 3. S-factor for the reaction ^Li(p,a)^He as a function of the in- 
cident cm. energy. The experimental points are from ((4SlfMarion56) . 
from (|4^)(Gemeinhardt66), from (|53)(Spinka71), from (jjll ) (Elwin79) . 
from (l5i)(Shinozuka79), from ([si) (Kwon89) and from @) (Engstler92) . In 
the top panel the solid curve represents our fit by polynomial expression with the 
adiabatic enhancement. The dashed curve corresponds to the bare S-factor. In the 
bottom panel the solid curve represents our fit by using Eq. ([5|), instead of the 
polynomial expression, and with treating the screening potential as a parameter. 

3.4 D(p,7)3He 



The S'-factor data of the reaction D(p,7)'^He from several measurements are 
shown with error bars in Fig. |H We performed the fitting procedure of the data 
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Table 6 

Fitting parameters and the 5*^(0) of the reaction ^Li(p, a)^He by using fitting pro- 
cedure Eq. ([5]). The first row is obtained by using the screening potential in the 
adiabatic limit, the second row corresponds to the zero screening potential, the last 
row is obtained by treating Ue as a fitting parameter. 



no/(o) 


do 


Ue (eV) 




56(0) (MeVb) 


0.11 


1.2 


175 


2.5 


3.3 


0.11 


1.2 





3.8 


3.5 


0.12 


1.1 


466 ± 31 


1.4 


3.1 



Table 7 

Fitting parameters of the reaction D(p, 7)^He. The three rows are obtained by us- 
ing a quadratic polynomial and from (l)(the first row), without enhancement (the 
second row) and with the adiabatic approximation (the third row). 



Sm (eVb) 


Si (b) 


52 (eV-^b) 


Ue (eV) 




0.20±0.07 


5.60±2.00 


3.10±1.10 


0.0 


(NACRE) 


0.261±0.006 


1.3±0.2 


12.0±1.0 


0.0 


2.7 


0.256±0.006 


1.4±0.2 


ll.Sibl.O 


40.8,0.0 


3.9 



Table 8 

Fitting parameters and the 5f,(0) of the reaction D(p, 7)^He by using the fitting 
procedure Eq. ([5]). The first and the second rows are obtained by the fitting with- 
out enhancement, and by using the screening potential Ue in the adiabatic limit, 
respectively. 



no/(o) 


ni/(0) 


do 


di 


Ue 


(eV) 


X// 


5fe(0) (eVb) 


2.5x10-^ 


3.4x10-*^ 


3.4 


4.0 




0.0 


2.7 


0.25 


2.3x10-*^ 


3.6x10"* 


3.7 


3.9 


40.8 


, 0.0 


2.8 


0.25 



in the incident energy region from 0.0025 (MeV) to 1 (MeV) using a quadratic 
polynomial without enhancement factor. The obtained fitting parameters are 
shown in the second row of Table [71 In (jlf) the same polynomial degree has 
been used but the low-energy data by (|.54|) were not available at that time. 



At the minimum incident energy in Ref. (1541 ) the screening enhancement is 
estimated to be 7% at utmost. This enhancement within the adiabatic limit 
is, again, estimated by using a linear combination of the even and odd states 
of the electronic wave function, reflecting the charge symmetry of the system 
as Eq. (USD where U^^^^^ and [/J^^)2 are replaced by 40.8 eV and 0.0 eV, 
respectively. 

The fitting parameters obtained using a quadratic polynomial with the adi- 
abatic enhancement are shown in the third row in Table [71 and it is shown 
with the dashed curve in the top panel in Fig. [H Because the enhancement 
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Fig. 4. 5-factor for the reaction D(p,7)'^He as a function of the incident cm. energy. 
The ex per imental points are from Issl) (Griffiths62) . from (^) ( Warren63) , 
from Isj) (Berman64) , from (|5^)(Wolfii67). from (j5^ )(Schmid95) and 
from (j54l )(Casena02). In the top panel the solid curve represents our fit by 
polynomial expression with the adiabatic enhancement. The dashed curve corre- 
sponds to the bare S'-factor. In the bottom panel the curves are same with ones in 
the top panel, but using Eq. ([5]) instead of the polynomial expression. 



is less than 7% even at the lowest measured incident energy, it changes in- 
significantly the zero-energy S'-factor: S'(0)=0.261 ± 0.006 (eVb) obtained by 
neglecting the enhancement differs only slightly from the bare S'-factor at zero- 
energy S'b(0)= 0.256 ± 0.006 (eVb) from our fitting procedure. S'b(0)= 0.256 ± 
0.006 (eVb) is slightly higher than the result from the i?-matrix fit SbiO)= 0.223±0.010 (eVb) 



17 



in Ref. (1211 ). which is obtained as a sum of Ml and El contributions. Limiting 
the fitting range from 0.0025 MeV to 1 MeV, Sb{0) is insensitive to the choice 
of the degree of the polynomial. However the Sb{0) obtained using polynomials 
is quite sensitive to the choice of both the upper and the lower limits of the 
fitting range. 

We performed the fitting procedures using Eq. ([5]). This fitting procedure with- 
out enhancement and with the adiabatic enhancement factor lead the fitting 
parameters and the S'b(O) in Table [H We have used the / = and 1 states. 
The obtained Sf,{0) are essentially the same for two cases and are in agreement 
with the extrapolations using quadratic polynomials. The radial parameters 
do and di for both cases are larger than one. This can be interpreted because 
the effective radius of deuteron is larger than the one given by the empiri- 
cal formula. The curve obtained by this fitting procedure for the adiabatic 
enhancement is shown in the bottom panel in Fig. HI 



3.5 3He(a,7)^Be 



The ^-factor of the reaction ^He(a,7)^Be has been investigated recently both 



by the activation (l32l : l33l ) and the prompt methods (j3J). The latter confirmed 
that there is no discrepancy between the obtained S'-factors by two different 
methods. They have discussed the electron screening enhancement factor in 
the adiabatic limit in (jsil ). but the S'-factor data has not been corrected by 
the effect. The S'-factor of the reaction ^He(a,7)''Be from several measure- 
ments are shown with error bars in Fig. O The fitting parameters in (jll) are 
shown in the first row of Table M The screening potential for the reaction 
^He(a,7)^Be is estimated in the same way with the reaction ^He(^He,2p)^He. 
The estimated enhancement at the minimum incident energy within the adi- 
abatic limit is 2% at utmost. We performed the fit of the data in the incident 
energy region from 0.1072 MeV to 1 MeV using quadratic polynomials with 
the adiabatic enhancement factor. The obtained fitting parameters are shown 
in the second row of Table M We have performed the same fit but without en- 
hancement factor and obtained the same Sb{0) as the one with the adiabatic 
enhancement. The obtained Sb{0) coincides with the one in ^ within the 
error. However Sb{0) is rather sensitive to the choice of the fitting range, espe- 
cially to the choice of the lower limit. Sb{0) is insensitive to the choice of the 
degree of polynomial in the selected fitting range. In the top panel in Fig. Owe 
have shown the results of the fitting using the adiabatic enhancement factor. 
The S'-factor at zero-energy Sfc(0)=0.49 ± 0.01 (keVb) from our procedure is 
in agreement with the result from the i?-matrix fitting 0.51±0.04 (keVb) in 



Ref. (1211). 



The fitting parameters from the fitting procedures with Eq. are shown in 
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Fig. 5. S'-factor for the reaction ^He(a,7)''Be as a function of the incident cm. 
energy. The experimental points are from (j6Q)(Parker63), from (|6ll ) (Kraewinkel82) . 
from (l6p ) fOsborne82) . from Ig^) (Hilgemeier88) . from (js^ ) (BemmererOG) and 
from (1341) (ConfortolaOT) . In the top panel the solid curve represents our fit by poly- 
nomial expression with the adiabatic enhancement. The dashed curve corresponds 
to the bare S-factor. In the bottom panel the curves are same with ones in the top 
panel, but using Eq. ([5]) instead of the polynomial expression. 



Table fTOl We have used the / =0, 2 partial-wave contributions. This procedure 
using the enhancement factors with U^^^^ = 246.8, 0.0 (eV) gives zero energy 
S'-factors 0.50 and 0.51 (keVb), respectively. Both are in agreement with the 
result from the polynomial fitting procedure. The radial parameters do and d2 
for both cases are larger than one. The curve obtained by this fitting procedure 
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Table 9 

Fitting parameters of the reaction '^He(a, 7)^Be. The two rows are obtained by using 
polynomial expression and from (']]) (the first row), with the adiabatic approxima- 



tion (the second 


row). 








St{0) (keVb) 


Si (b) 


52 (keV-^b) 


Ue (eV) 




0.54 ± 0.09 


-0.52 


-0.52 


0.0 


(NACRE) 


0.50 ± 0.01 


-0.35 ± 0.04 


0.13 ± 0.03 


246.8 


0.0014 



Table 10 

Fitting parameters and the Sb{0) of the reaction ^He(Q,7)''Be by using the fitting 
procedure Eq. ([5]). The first and the second rows are obtained by the fitting with- 
out enhancement, and by using the screening potential Ue in the adiabatic limit, 
respectively. 



no/(o) 


n2/(o) 


do 


d2 


Ue (eV) 




Sb{0) (keVb) 


5.1x10""^ 


4.9x10^^^ 


4.0 


2.5 


0.0 


2.3 


0.51 


5.3x10-'^ 


4.9x10"'^ 


4.0 


2.5 


246.8 


2.3 


0.50 



with the adiabatic enhancement is shown in the bottom panel in Fig. [51 The 
obtained zero energy 5'-factor is smaller than 5'(0) = 0.560 ± 0.017 (keV) 
in (I34J). This is because their result is obtained by a normalization to their 
data. 



3.6 ^Be(p,7)8B 



The reaction ''Be(p,7)^B is a key process to produce the high energy solar 
neutrinos through the /9-decay of ^B. The S-factor of this reaction is inves- 
tigated intensively by many groups by means of the direct capture (DC) re- 
action (1351 : 1361) ■ the indirect Coulomb dissociation (CD) method (1371 : l39l) and 
the asymptotic normalization coefficients (ANCs) (1641 ). It was claimed that the 
experimental data of the S'-factor by CD experiments gives a steeper energy 
dependence than that by DC experiments in the low-energy region and the 
lower zero-energy S'-factor as an average (1361 ) . Recently this inconsistency has 
been resolved by reanalyzing data by the Coulomb dissociation method (39^. 



For the purpose of the extrapolation of zero-energy S'-factor, both the DC and 
the CD experiments above mentioned use the microscopic cluster model (IHsl ) 
and give the zero-energy S'-factor 22.1 ± 0.6(expt.) ± 0.6(theor.) (eVb) [21. 4± 
0.5(expt.) ± 0.6(theor.) (eVb) as a mean of all modern direct measurements] 
and 20.6 ±0.8 (stat.) ± 1.2 (syst.) (eVb), respectively. Let us remind you that 
the purpose of this paper is not a precise determination of the S'-factor but to 
see the effect of the electron screening on the determination of the S'-factor. 
So that we rather use the consistent approaches with the other reactions than 
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Fig. 6. S'-factor for the reaction '^Jie(p,y)^B as a function of the incident cm. energy. 
The experimental points are from (|36l ) (JunghansOS) and from (,39i) (Schuemann06) . 
In the top panel the solid curve represents our fit by polynomial expression with 
the adiabatic enhancement. The dashed curve corresponds to the bare S'-factor. In 
the bottom panel the curves are same with ones in the top panel, but using Eq. ([5|) 
instead of the polynomial expression. 



employ a special treatment for this reaction. However we, at least, need to 
include the resonances to analyzing the DC data (1361 ). For this purpose we use 
the Breit-Wigner formula Eq. (Q. Moreover it is well known that this reaction 
has a low-energy bound state in the (1661 : 1671 ). To take into account this 
state, we use 
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Table 11 

Fitting parameters and Sb{0) of the reaction ^Be(p,7)^B by using the polynomial 
Eq. ()16p . The four rows are obtained from (J,j)(the first row), without enhance- 
ment (the second row) , with the adiabatic approximation (the third row) , and with 
the screening potential Ue= 900 (eV) (oS) (the fourth row). 



S-i (eV^b) 


So (eVb) 


Si (b) 


S2 


Ue (eV) 




sm 




21. ± 2. 


18. 


38. 


0.0 


(NACRE) 


21. ± 2. 


0.6 ± 0.2 


15.7 ± 0.6 


8.3 ± 0.4 




0.0 


0.3 


20.8 


0.7 ± 0.2 


15.7 ± 0.6 


8.3 ± 0.4 




222.0 


0.3 


20.5 


0.7 ± 0.2 


16.0 ± 0.6 


8.1 ± 0.4 




900.0 


0.3 


19.7 



Table 12 

Fitting parameters and the ^^(O) of the reaction ''Be(p, 7)^6 by using the fitting 
procedure Eq. ([5]) with resonant terms. The three rows are obtained by the fitting 
without enhancement (the first), and by using the screening potential Ue in the 
adiabatic limit (the second), and with the screening potential Ue= 900 (eV) (the 



third), respectively. 


no/(o) n2/(o) 


do 


d2 


Ue (eV) 




^0 (eVb) 


2.2x10"^ 2.5x10-^ 


0.9 


0.2 


0.0 


0.4 


22.4 


2.2x10"*^ 2.3x10"^ 


0.9 


0.2 


222. 


0.4 


22.3 


2.2x10-*^ 2.0x10"^ 


0.9 


0.2 


900. 


0.5 


22.1 



S,{E) = + So + SiE, (16) 



where Eb = 0.1375 (MeV) (1661 ) . in place of the polynomial expression. We 



make a special mention of ^Be metallic target being used in the experi 



ment (1361 ) . The experimental data of the 5'-factor in Ref. (1361 ) is fitted with 
the Breit-Wigner single-level resonance formula for l"*" and 3"*" resonances plus 
Eq. (fT6!l . We use the resonance parameters in Ref. (j36l ) and obtained 5';,(0) = 
20.8 (eVb) with xt =0.3 from the fitting procedure without enhancement fac- 
tor. Assuming the adiabatic enhancement, the screening enhancement factor is 



of the order of 1% at the minimum incident energy of the experiment in (1361 ). 
By making use of the enhancement factor with the adiabatic screening po- 
tential f/i^^) = 222.0 (eV), we obtain ^^(0) = 20.5 (eVb) with xl =0.3. The 
corresponding bare S-factor is shown with the dashed curve in the top panel in 
Fig. [61 The difference between two Sb{0) using different screening potentials is 
less than 2%. Considering the isotopic independence of the electron screening, 
we, tentatively, use the screening potential obtained by the measurement of 
the reaction ^Be(p, a)^Li: f/e=900±50 (eV) (jisi ). The fitting procedure gives 
Sb{0) = 19.7 (eVb) with xl = 0.3. This is 4% smaller than the former two 
results. We summarize the fitting parameters in above fitting procedures in 
Table M 
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We perform the fitting procedures of the experimental data in Ref. (1361 ) using 
Eq. ([5]) including l"*" and 3"'" resonances plus another resonance with negative 
energy -0.1375 (MeV), which corresponds to the pole term in Eq. fll6p . 



c,ie^»no^(E) c,2e^°no;(E) Crse'^myjE) 

{E + 0.1375)2 (E- 0.630)2 ^y^A (E- 2.183)2 + vyi ^ 

where 0^1,0^2 and 0^3 are scaling factors in Eq. (fTOj) . Our fitting procedures 
using different screening potentials give common Cri = 5. x 10"'' (MeV2), 
Cr2 = 2. X 10-9(MeV2) and = 2. x 10"^ (MeV2). In the bottom panel 
of Fig. E] the thin dotted and dot-dashed curves show the contributions from 
resonances, including the negative-energy resonance, and the non-resonant 
part, respectively. We have used the I = and 2 partial- waves. 

The fitting procedure using Eq. (|T7j) without enhancement, with the adiabatic 
enhancement factor, and with Ue= 900 (eV) lead the fitting parameters and the 
Sb{0) in Table [T2l Despite of the difference of the utilized screening potentials, 
neither the obtained zero-energy S'-factor nor xl does have much differences. 
The values 112/(0) are much larger than no/(0), i.e., the rf-wave contribution is 
dominant in the energy region investigated. The radial parameter do is smaller 
than one, and di is smaller than do for all three cases. The curve obtained by 
this fitting procedure for the adiabatic enhancement is shown in the bottom 
panel in Fig. [61 The obtained Sb{0) using different screening potentials are in 
accordance with the result with the microscopic cluster model ^) S';,(0)=22.1 
± 0.6(expt.) ± 0.6(theor.) (eVb) in Ref. (1361 ) within the error-bar. 

From the tree results of all the radiative capture reactions considered, D(p,7)^He, 
^He(a, 7)''Be, and ''Be(p, 7)^B, one can say that the obtained Sb{0) is insen- 
sitive if the screening enhancement, within the adiabatic approximation, is 
taken into account or not. For all the reactions investigated in this paper 
the fitting procedure with the polynomial expression is more sensitive to the 
difference of the screening potential than the fitting procedure with Eq. ([5]). 



4 Conclusions 



We discussed the bound electron screening corrections to the bare S'-factors 
of the reactions in PP-chains. Our approach is based on fitting procedures 
of the experimental data. For this purpose we employed two different fitting 
procedures: one is the conventional polynomial expressions and the other in- 
cludes explicitly the contribution of the nuclear interaction and based on the 
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statistical model to describe exit channels. The later fitting procedure works 
especially well for the reactions that have a dominant s-wave entrance chan- 
nel component. We have applied different types of screening enhancements: 
in the adiabatic limit, determined through a fit and larger than adiabatic 
screening potentials, as well. Prom the tree results of all the transfer reactions 
considered, ^He(^He,2p)^He, ^Li(p, Q;)^He, and ^Li(p,a)^He, the enhancement 
by the screening is crucial, in the sense that the fitting procedure without 
enhancement gives xt larger than the others in which the screening enhance- 
ment is taken into account. However the obtained Si,{Q) is insensitive to the 
magnitude of the screening potential. Especially for the radiative capture reac- 
tions D(p,7)^He and ^Be{p, 7)^B, the screening correction within the adiabatic 
approximation has been considered for the first time. However the results sug- 
gest that the obtained 5'{,(0) is insensitive whether the screening enhancement, 
within the adiabatic approximation, is taken into account or not. Making a 
comparative study of the bare (^-factors obtained by two- ways of fitting proce- 
dures using different screening enhancement factors, we found that all Sb{E) 
coincide within the standard errors. Sb{E) is, practically, insensitive to the 
magnitude of the screening potential. 

The authors acknowledge Prof. S. Kubono for the suggestion of the problem 
and valuable comments. One of us (S. K.) thanks Dr. H. Costantini and Dr. 
R. G. Pizzone for stimulating discussions and for providing us experimental 
data. 
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